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ABSTRACT 

We construct a theoretical model for low-temperature crystallization of amor- 
phous silicate grains induced by exothermic chemical reactions. As a first step, 
the model is applied to the annealing experiments, in which the samples are 
(1) amorphous silicate grains and (2) amorphous silicate grains covered with an 
amorphous carbon layer. We derive the activation energies of crystallization for 
amorphous silicate and amorphous carbon from the analysis of the experiments. 
Furthermore, we apply the model to the experiment of low-temperature crys- 
tallization of amorphous silicate core covered with an amorphous carbon layer 
containing reactive molecules. We clarify the conditions of low-temperature crys- 
tallization due to exothermic chemical reactions. Next, we formulate the crystal- 
lization conditions so as to be applicable to astrophysical environments. We show 
that the present crystallization mechanism is characterized by two quantities: the 
stored energy density Q in a grain and the duration of the chemical reactions r. 
The crystallization conditions are given by Q > Q m m and r < r coo i regardless of 
details of the reactions and grain structure, where r coo i is the cooling timescale 
of the grains heated by exothermic reactions, and Qmm is minimum stored en- 
ergy density determined by the activation energy of crystallization. Our results 
suggest that silicate crystallization occurs in wider astrophysical conditions than 
hitherto considered. 

Subject headings: planetary systems: protoplanetary disks- comets: general- meteors, 
meteoroids 
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Introduction 



Infrared ob servations of various comets suggest the existence of crystalline silicates 



in their nuclei (JBregman et al. 



19871 : 



Molster et al. 



19991 ) . In contrast, silicate in the 



inters tellar medium is almost or entirely amorphous, namely, free from crystallites (ILi et al. 



20071 ). It is, therefore, popular to assume that crystalline silicates observed in comets 
formed in the inner solar nebula (e.g., Gail 2001; Bockelee-Morvan et al. 2002; Harker and 
Desch 2002; Keller and Gail 2004). This is out of harmony with the composition of the gas 
in cometary comae , observations of which indicate the p reservation of int erstellar i c es in the 



cold outer nebula (JBiermann et al 



1982 : 



Mumma 



19961 ). As noticed by 



Huebner! ( 120021 ) 



none of the thermal mechanisms proposed for silicate crystallization allows comets to retain 
the interstellar composition of ices in their nuclei. 

Besides comets, infrared spectra tell us the presence of crystalline silicates in 
various kinds of object such as AGB stars, post-AGB stars, red supergiants, Herbig 
Ae/Be stars, and protoplanetary disks around young stellar objects. It is known that 
crystallization of amorphous silicate due to annealing requires te mperature T above 1000 K 



f lHallenbeck fc Nuth 



1998 : 



Fabian et al. 



2000 : 



Murata et al. 



20071 ) . It is generally considered 



that silicate does not c rystallize be 



forsterite composition fjSpeck et al- 



ow it s glass transition temperature, which is 990 K for 



20081 ) . However, infrared spectra of dust shells around 
evolved oxygen-rich stars exhibit the presence of several emission features of crystalline 
silicates in their cool dust shells (T < 300 K), while sili cates condense in d usty outflows 



with amorphous structure as evidenced by the spectra (IWaters et al 



19961 ). Crystalline 



silicates are also observed in ULIRGs (ultr aluminous infrared galaxies), implying that 



they are located in the cool, outer regions ( iSpoon et al.l 120061 ) . From these observational 



results, several researchers clai m a necessity of a yet unknown crystal 



low-temperature environments (IWaters et al. 



1996 



Molster et al. 



1999 



i zation process in 



2001 



Spoon et al. 
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20061 : IWatson et all 120091 ). 



Physically, crystallization is a re-arrangement of atoms irregularly placed in a solid 
so that the atoms occupy lattice sites; in such a configuration the solid attains the lowest 
internal energy. However, the atoms must overcome the energy barrier, namely, activation 
energy of crystallization E. In annealing, a fluctuation of thermal energy is used to 
overcome the energy barrier. However, a flu c tuatio n of the thermal energy is not necessarily 



a requisite for crystallization. 



Carrez et al. 



Kimura et al. 



(I2008H ) showed that 



(120021 ) and 

amorphous silicate with forsterite composition crystallizes at room temperature when they 
are irradiated by electrons in a transmission electron microscope (TEM). 

Recently, we have proposed another mechanism of low-temperature crystallization due 
to exothermic chemical reactions of reac tive molecules in an organ i c refractory mantles 



surrounding an amorphous s ilicate core (JYamamoto fc Chigai 



2009 



Yamamoto et al 



2005 



Yamamoto fc Tanaka 



2009J ). Once the reactions are triggered by moderate heating above 



a few hundred kelvins, chain reactions of reactive molecules contained in the organic mantle 
release heat that crystallizes the surface layer of the amorphous silicate core. We estimated 
the degree of crystallinity to be 0.4 to 20% in volume for the concentration of reactive 
molecules of 1 to 10 %. It was shown that the degree of crystallinity was sufficient to 



reproduce the ob served strengt 



i of in 



comet ary comae (JKimura et al 



2008a 



rared features characteristic of crystalline silicates in 



2009|) 



and 



Nevertheless, we notice that the previous model of 



Yamamoto et al. 



Yamamoto fc Chigail (120051 ) 



(120091 ) is incomplete, because the model permits crystallization at 
temperatures above the melting point, although, in fact, the existence of crystalline silicate 
is not allowed in these temperatures. If one excludes crystallization at temperatures 
higher than the melting point, one would have a lower volume fraction of crystalline 
silicate. Furthermore, their model neglects the effect of finite particle size, namely, neglects 
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accumulation of the heat of chemical reactions in the particle of a finite size. In consequence, 
we expect that the particle-size effect increases the degree of crystallinity. These two 
points should be addressed in the model to properly examine the plausibility of nonthermal 
crystallization due to exothermic chemical reactions. 

In this study, we construct a model for low-temperature crystallization of an amorphous 
silicate core coated with a layer of carbonaceous material, by taking into account the 
melting point and a finite particle size (in section 2). In view of the analyses of cometary 
dust and primitive interplanetary dust particles, we consider a grain having an amorphous 
silicate core and a mantle of carbonaceous material. Amorphous carbon is one of the main 



components in t he carbonaceous material. (IKeller et al. 



Jessberger 



1996 



2000 



Kissel fc Krueger 



1987 



19991 ). We determine the activation energies of crystallization in amorphous 



silicate and amorphous carbon using the revised model in comparison with two previous 
crystallization experiments, in which particles of (1) amorphous silicate and (2) amorphous 
silicate covered with a carbonaceous layer are annealed. These evaluations enable us to 
derive quantitative crystallization conditions. We demonstrate the validity of the model 
by its applicatio n to crystalli z ation experiments due to exothermic chemical reactions 



demonstrated by 



Kaito et al. 



( 120071 ) in section 3. Finally we formulate the crystallization 
conditions applicable to various astrophysical environments in section 4. We discuss the 
feasibility of nonthermal crystallization in astrophysical environments and summarize our 



conclusions in section 5. 
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Model 



2.1. Basic equations 



Cometary dust and interstellar dust are well modeled as aggreg ates of small grains 



having an amorphous silicate core and an organic refractory mantle ( IKimura et al 



20031 ) 



In-situ measurements of mass spectra for dust in a coma of comet lP/Halley show evidence 
of such a core-mantle structure, in which a core and a mantle ar e composed of silicate and 



organic refractory material, respectively (IKissel fc Krueger 



19871 ). A cluster of submicron 



grains with amorphous silicate enclosed within carbonaceous material such as amorphous 
carbon rich in o rganic material is common in primitive interplanetary dust part icles of 



cometary origin (IKeller et al 



1996 . 



2000 : 



Kissel fc Krueger 



1987 



Jessberger 



1999). Because 



silicate crystallization in each constituent grain of the cluster occurs independently, we 
hereafter consider a single particle having a silicate core covered by amorphous carbonaceous 
material. We expect that reactive molecules have been formed by the exposure of ultraviolet 
radiations and high energy particles in molecular clouds and stored in the organic layer of 
interstellar dust. 

Once the dust is heated up, the diffusion rate of reactive molecules increases and 
chemical reactions are triggered in the mantle layer. The energy released by the reactions 
raises the temperature and in turn expedites further reactions. As a consequence of the heat 
flow, the amorphous silicate core crystallizes. If the mantle layer is composed of amorphous 
carbon, the amorphous carbon may also be graphitized (Kaito et al. 2006, see also Speck et 
al. 2009). We take into account two types of crystallization of the silicate core. 

The timescale of heat transfer in a particle is estimated to be 



Theat = r 2 / X = 3 X 10" 



X 



v. 



10 3 cm 2 s x 
for the particle radius of r = 50 nm. The values of x range from 4 to 20 x 10~ 3 cm 2 s _1 de- 
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pend i ng on the temperature, 



2009 



Hofmeister et al. 



compositions, and degree of crystallization ( IWhittington et al. 



20091 ). Thus we find 7"h ea t < 10~ 8 s, which is much shorter than 
the crystallization timescale > 10~ 5 s (see Appendix A). Consequently, it is safe to assume 
that a particle is isothermal during the crystallization process. 

The energy budget of the particle at temperature T is described by 

A ra d — A co u, (2) 



4 dT A 

-7r(o + hfp A c p — = -tt[(o + hf - a 3 }e + T si + T c 



where a is the radius of the silicate core, h is the thickness of the carbonaceous mantle, 



Pd(— 3.0 gem 3 ) is the mean density of the particle, a nd c v {= 1.5 x 10 



1995 



ergg 



Lide 



1 K~ 1 ) is the 



19961). The 



mean specific heat of amorphous carbon and silicate (jNavrotskyi 
first term on the right-hand side in Eq. (T2J) is a heating term due to exothermic chemical 
reactions and r si is a heating rate due to deposition of the latent heat of crystallization of 
amorphous silicate, and T c is that of amorphous carbon. For cooling processes, we take 
into account two cooling mechanisms due to thermal radiation A ra d and collisions of gas 
molecules surrounding the particle, A co n. We assume a first order reaction, namely, the 
main reactions are those of reactive molecules contained in the mantle with the mantle 
materials. 



Then the time variation of the number density n^ of the reactive molecules is expressed 



by 



driA nx 

~dT ~ T ; 



(3) 



where r is a timescale of the reactions. The heating rate e by the reactions is given by 



drix kQ ( t 

e = — — q x = — exp — 

at t \ r 



(4) 



Here, Q is the energy density stored in the carbonaceous layer given by Q = nj^^/k, 
where ?t,a,o is the initial number density of reactive molecules, and q T is energy release per 



reaction. The cooling rate due to thermal radiation, A ra( j, is given by 

A rad = 47r(a + / i )V B e(T 4 -T 4 ), (5) 

where e is the efficiency of thermal emission from the surface of the carbonaceous mantle, 
T is the ambient radiation temperature, and <jb is the Stefan-Boltzmann constant. The 
cooling rate A co n due to collisions of ambient gas molecules of temperature T is given by 

A co ii = 2n{a + h) 2 nvk{T-T ). (6) 

Here, n is the number density of the gas molecules, k is the Boltzmann constant, and 



v = \J 8/cT //im a is their mean thermal velocity, where // is the mean molecular weight and 
m a = 1.6 x 10~ 24 g is atomic mass unit. 

We assume that the crystallization proceeds from the interface between the core and 
the mantle as is supported by crystallization experiments of grains having a silicate core 
and an amorphous carbon mantle (see section 2.2). This means that a crystal growth 
occurs through heterogeneous nucleation on the interface. Namely, crystallization of an 
amorphous silicate core proceeds inward in the core from the interface and graphitization 
of an amorphous carbon mantle proceeds outward in the mantle from the interface. 

Denoting the distance of the silicate crystallization front by a s ; and the graphitization 
front from the center of the particle by o c , equations of the crystal growths in the silicate 
core and the carbonaceous mantle are given, respectively, by 

<?i,siAT £ 



da s \ 1/3 / E R 

-il! ^ si exp -- 



dt S1 S1 "V kT 



1 — exp 
- exp 



kT 2 
<?i, c AT c 



(7) 



kT 2 

where E^ (i = silicate (si), carbon (c)) is activation energy of crystallization, f2j is volume 
of the unit cell, Vi ~ 10 13 s _1 is vibration frequency of the lattice. In the square bracket, 
gi S i is latent heat of crystallization from the melt per unit cell of silicate, q\ yC is that of 
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graphitization from amorphous carbon, ATj = T c .j — T is the supercooling, in other words, 
the temperature difference between a temperature T ei of the transition from the melt or 
the amorphous state to the crystal in equilibrium and the temperature T in concern. The 
the square-bracket factors indicate the correction near the melting temperature and almost 
equal unity at temperatures except near the melting temperature. The heating rates due to 
deposition of the latent heat of crystallization are expressed by 



-4vra 2 gi ' si dasi r ■ - «~» 2Ql ' cdac 



fL; dt 



A-nai 



Q r dt 



(9) 



2.2. Determination of activation energy of graphitization 



First we shall derive the activation energies of crystallization in amorphous silicate and 
amorphous carbon, (i.e., E si and E c ), by applying the model described in the previous section 
to crystallization experiments on submicrometer-sized grains having an amorphous silicate 
core and an amorphous carbon mantle. The values of the activation en ergies will be used in 



KaitoetaL 



(12007]). 



the analysis of the more sophisticated experiment of crystallization by 

The upper panel of F ig. [T] shows a schematic p i cture of the crysta llization experiments of 



amorphous silicate by 



Kaito et al. 



(12006 . 



20071 ). 



Kaito et al. 



(120061 ) performed an annealing 



experiment on amorphous silicate particles (~ 50-100 nm in diameter) covered with an 
amorphous carbon layer (~ 20-30 nm). In this experiment (which we call Exp. 1 hereafter), 
the crystallization of the amorphous silicate core of forsterite composition was observed at 
the ambient te mperature of 870 K in vacuum, which is 2 00 K lower than the case for bare 



silicate grains ( IFabian et a 



layer was also graphitized. 



2000 



Kamitsuji et al. 



Kaito et al. 



20051 ). The amorphous carbonaceous 



(120061 ) also found that the crystallization in both the 



silicate core and the carbon mantle proceeded from the interface between the core and the 
mantle. It might seem curious that the "crystallization temperature" of 870 K is lower than 
the glass transition temperature ~ 990 K for amorphous silicate of forsterite composition 
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( jSpeck et al.ll2008l ). Actually, crystallization of the amorphous silicate core occurs above 
the glass transition temperature in Exp. 1 as well. Indeed, our analysis of Exp. 1 shows 
that crystallization of amorphous silicate occurs around 2000 K in a short time because of 
a rapid temperature elevation in the grain due to graphitization of the amorphous carbon 
mantle followed by the cooling (see Fig. [3]). 

Sublimation of the grains is negligible even though they suffer temperatures as high 
as 2000 K. We shall show this by comparing sublimation and cooling timescales. The 
sublimation timescale is estimated by 

1 \ fd(a + h) \ ] -1 _ p d (a + h) 



\a + hj V dt 
where F e is the sublimation rate given by 

F e = aP sat (T) 



m 



11) 



2vrfcT 

with a being the evaporation coefficient. The vapor pressure P sat is approximated by the 
Clapeyron-Clausius relation expressed by P sat = Poexp{—H/(kT)}, where Po is a constant 
H is latent heat of sublimation, and m is mass of the subliming molecules. On the other 
hand, the cooling timescale including both collisional and radiative coolings is given by 

An (a + h) 3 p d c p T 2p d c p (a + h) 



7~cool 



3(Acou + A rad ) 3(nkv + 2a B £ e T c 3 ) 



(12) 



for T > Tn. 



Figure [2] compares the sublimation and cooling timescales of a grain coated with a 
carbon mantle; the cooling timescales are estimated for the ambient gas pressure of 1 atm 



and vacuum. For the set up of the Exp. 1 (jKaito et al 



and p = 29 (jSatoh 



20061 ). we put a = 50 nm, h = 20 nm, 



l2002l ). Because of the core- mantle structure of the grain, sublimation of 
the ca rbon mant le would occur first, for which H/k = 9.5 x 10 4 K, P = 1.3 x 10 15 dyn 



cm 2 (ILide 



1996|), an d oc = 1. We calculated the emissivity to be e = 0.01 for a core-mantle 
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grain of a = 50 nm and h = 20 nm using Mie theory. Figure [2] indicates that the sublimation 
timescale r su b is longer than the cooling timescale r coo i for T < 3700 K in vacuum and 
T < 5200 K at 1 atm. The temperature of the grain experienced during the heating stage 
is substantially lower than these temperatures (see Fig. [3] for example), thus sublimation is 
negligible. This is consistent with the results of the experiments that the grains suffer little 
sublimation. 

In this section, we apply the model to Exp. 1 to derive the activation energy of 
amorphous carbon. We solve the energy equation of the particle given by Eq. ([2]) without 
chemical reactions in the carbon mantle (e = 0) and cooling due to collisions of surrounding 
molecules since Exp. 1 was performed in vacuum (A co ii = 0). 

There are two cases of crystallization of the silicate core in the present mechanism, 
namely, crystallization from amorphous solid or the melt (liquid). The latter case occurs if 
the temperature of the particle exceeds the melting point by the large heat deposition as 
will be shown later. In this case, the melt crystallizes in a subsequent cooling stage when 
the temperature becomes lower than the melting point. In general, the activation energy for 
crystallization E si has different values for crystallization from amorphous solid and liquid. 
Thus we set 



E s 



E K 



E s ; s for t < t m , 



E s u for t > t r 



(13) 
(14) 



where t m is time when the particle temperature reaches the melting point, E 8 i >8 and E s ^\ are 
the activation energies for the amorphous silicate and that of the melt, respectively. We take 



E S i,s/k = 39000 K, which was obtainec 



bare amorphous silicate (JFabian et al 



from the analysis of crysta llization experiments of 



2000 



Kamitsuji et al 



20051 1 (see also Appendix [A} . 



On the o ther hand, the valu e of the activation energy of liquid E si ^/k is taken to be 



23000 K (ITanaka et al. 



20081 ). For other parameters, we adopt q\ jS i/k = 13800 K ( jNavrotsky 
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19961 1. fi si = 7.3 x 10~ 23 cm- 3 (INavrotskv 



Vr = U s : 



1.2 x 10" s" 1 flFabian et al 



1995 ). qijk = 5100 K (IKaitol 120070. T oai = 2160 K (INavrotskvlll995h . T P . r = 3915 K (ILide 



2000). 



19951), ^c = 8.0 x 10- 24 cm- 3 (ILidelll996h . and 



Figure [3] shows a typical feature of temporal variations in the temperature T of the 
particle and the thicknesses / s ; = a — a S [ of the crystallized layer in the silicate core, and 
that of the graphitized layer, l c = a c — a, in the carbon mantle at To = 870 K. Here we take 
E c /k = 23000 K, because this value reproduces the experimental results as will be shown 
later. 

The behavior of the crystallization process is described as follows. First, graphitization 
proceeds in the carbon mantle, starting from the interface between the mantle and the 
silicate core. During the growth of the graphite layer, the latent heat deposited is gradually 
accumulated in the particle because the timescale of radiative cooling ( > 1 s) is longer than 
the growth timescale of the graphite layer (~ 0.1 to 1 s) at T < 1000 K. This leads to a 
gradual increase in the temperature of the particle. The temperature rises by about 100 K 
from 870 K to 1000 K until the thickness l c of the graphite layer becomes a few nano-meters. 
A slight increase in the temperature at t = 0.076 s leads to rapid graphitization. For 
example, the growth rate of graphite at 1000 K is 20 times larger than that at 870 K ( 
i.e., 5 x 10 -5 cms -1 ). This results in a rapid deposit of the latent heat of graphitization and 
in turn the temperature rise higher than the melting point of forsterite. 

The rapid rise in the temperature is a sort of a positive feedback process, in which the 
graphitization in the carbon mantle releases the latent heat deposition and this heating 
accelerates graphitization. Once the graphitization of the mantle completes, the particle 
cools by thermal emission. When the temperature decreases below the melting point of 
silicate, crystallization of the silicate core begins and proceeds inward starting from the 
interface between the silicate core and the carbon mantle. The silicate crystallization stops 
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when the crystallization of amorphous silicate reaches the center of the silicate core. The 
timescale of silicate crystallization is about 10~ 3 s. After the complete crystallization of the 
amorphous silicate core, the particle cools further by thermal radiation. 

Figure H] shows crystallization features with varying the values of the activation energy 
E c of gr aphitization a t two ambient radiation temperatures of T = 770K and 870K. The 



result of 



Kaito et al. 



(120061 ) indicates that bare amorphous silicate does not crystallize at 
T = 770 K but does at T = 870 K. This fact constrains the value of E c . Figures H](a)-(c) 
show that the E c - value of E c /k < 25000 K explains crystallization of amorphous silicate 
at T = 870 K, while Figs. 13(d)- (e) show that E c /k > 21000 K yields no crystallization 
at T = 770 K. Thus, the activation energy of amorphous carbon must lie in the range of 
21000 K< E c /k < 25000 K. 

The range of the .E c -value is also estimated by comparing the growth timescale of 
graphite with the cooling timescale 7" coo i- If the growth time is shorter than the cooling 
time t coo i, amorphous carbon can be graphitized. Therefore, the condition of the graphite 
growth is: 

h 



da c /dt 
which leads with the use of Eq. (jSJ) to 



< Tcooi, (15) 



K < Tln l ^cooi^c 3 j = 23000K; (16) 

where we set T = 870 K and r coo i = 1 s in obtaining the value of E c on the right-hand side. In 
the same way, the condition of no graphite formation at 770 K is given by E c /k > 19000 K. 
The range of E c thus estimated is consistent with the range of E c constrained from Fig. HI 

From detailed calculations as shown in Fig. HI we find that the range of E c to explain 
the experimental results is 

21300 K< E c < 23600 K. (17) 
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In what follows, we use E c = 23000 K as a standard value of the activation energy of 
graphitization. 



3. Low- Temperature Crystallization in Laboratory Conditions 
3.1. Example of low-temperature crystallization 



Kaito et al.l ( 120071 ) performed another experiment on crystallization of amorphous 
silicate particles coated with amorphous carbonaceous layers (see the lower panel of Fig. [[]). 
In contrast with Exp. 1, however, the carbon coating was done in a CH 4 atmosphere of 
its pressure of 10~ 3 Torr (which we call Exp. 2 hereafter). The sample was observed by 
TEM after exposure to the air of the pressure of 1 atm. They found crystallization at room 
temperature (~ 300 K) near the interface of the silicate core and the carbon layer. The 
crystallization time is less than one minute in the air (Kaito, private communication). They 
consider that heat of oxidation of methane reacting with oxygen in the air graphitized the 
amorphous carbon layer and the latent heat of graphitization in turn induces crystallization 
of the amorphous silicate core of the particles. 

To clarify the conditions of crystallization due to exothermic reactions, we solve the 
equations of the model described in section 2 taking into account exothermic reactions and 
cooling due to collision of ambient air molecules at T = 300 K and the pressure of 1 atm. 
Figures and [T5] show a few examples, indicating time variations of the temperatures and 
the thicknesses of the crystalline silicate and graphite layers. In these figures, the radius of 
a silicate core and the carbon mantle thickness are the same as in Exp. 1, (i.e., a = 50 nm 
and h = 20 nm). Figure [5] are those for r = 5 x 10~ 8 s and (a) Q — 1.1 x 10 27 Kcm -3 , 
(b) 1.0 x 10 27 Kcm~ 3 , and (c) 0.9 x 10 27 Kcm~ 3 . In contrast to the case of Exp. 1, heat of 
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reactions of molecules (CH4 in Kaito et al. (2007)) contained in the carbon mantle with 
air leads to gradual rise in the temperature at the first stage and triggers graphitization of 
the carbon mantle. Heat released in graphitization leads further temperature rise above the 
melting point of silicate as seen in Figs. G^a)-(b). 

Silicate crystallization is observed in the subsequent cooling below the melting 
temperature. Note that the timescale of these processes is much shorter than that in 
Exp. 1 because the cooling of the particle is determined by the collisions of air molecules in 
Exp. 2. For Fig. E(a)-(b), the thickness of crystalline layer of silicate is about 7 nm. The 
thickness of the silicate crystalline layer is consistent with the experimental results, which 
are shown by the bars with arrows in the bottom panels in Fig. [5j On the other hand, no 
silicate crystallization is observed in Fig. [5]^c). This is because the stored energy density 
Q is insufficient for the temperature to rise above the melting temperature. From the 
calculations with varying Q by a finer step, we find that Q > 1.0 x 10 27 Kcm -3 is necessary 
for silicate to crystallize. In Fig. [61 we set Q = 1.0 x 10 27 Kcm -3 and varied the reaction 
timescale: (a) r = 5 x 10~ 9 s, (b) 5 x 10~ 8 s, and (c) 2 x 10~ 7 s. We find that r < 5 x 10~ 8 s 
is necessary for silicate crystallization. Both conditions of Q > 1 x 10 27 Kcm~ 3 and 
r < 5 x 10 -8 s are necessary for silicate crystallization. 

3.2. Analysis of Exp. 2 

We give analytical consideration of Exp. 2 and estimate a minimum value of the stored 
energy density Q for inducing crystallization of the silicate core in the conditions of Exp. 2. 
Since cooling of the particle is determined by collisions of air molecules in Exp. 2, the 
timescale of the cooling is given by 

47r(a + h) 3 p d CpT 2p d c p (a + h) 

r cooi = 5-r 5— rz (T > T ), (18) 

3A co u Znkv 



-16- 



which is about 10~ 6 s at T ~ 1000 K. This timescale is much shorter than the timescale 
of the crystallization of amorphous silicate, which is longer than 10~ 5 s (see Appendix |A~|) . 
Therefore, crystallization from the solid amorphous silicate is impossible in the conditions 
of Exp. 2. However, once amorphous silicate melts, crystallization is possible since 
crystallization from the melt is easier than from the amorphous solid. The timescale of 
crystallization from melting particles r cry is estimated from 



a 



T n 



(19) 



cry da S i/dt 

of which the minimum value is about 10 -6 s at T ~ 2000 K. This timescale is on the same 
order of magnitude as the cooling time in Exp. 2. 

As shown in the previous section, graphitization of amorphous carbon expedites 
crystallization of the silicate core. The amorphous carbon mantle crystallizes after the 
temperature increases due to chemical reactions in the mantle. From Eq. ([2j), the time 
variation of the temperature before graphitization (i.e., r si = T c = 0) is given by 



T = Tn + 



f^^c '< cool 



PdCp(r coo i - r) 



a + h 



-t/r 



-t/r 



(20) 



where r coo i is timescale of the collisional cooling given by Eq. (|T5|) . Here we ignore the 
radiative cooling because the cooling is mainly due to collisions of air molecules in Exp. 2. 
The temperature given by Eq. (|20|) increases with time and reaches a maximum at 

ln(r coo i/r) 



t 



1/t- 1/rcool' 



(21) 



The maximum temperature T c is given by 

iv\o£ Ti COOl 



T P -T„ 



PdC p (r cool - r) 
kQ 



1- 



a + h 



''"cool 



t/(t coo i— t) 



T 



'cool 



PdCp 



a + h 



(r < r C ool)- 



(22) 
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Note that T c is independent of r and r coo i and determined only by Q. Equation (122]) gives 

the minimum stored energy density Q m i n given a temperature T c that induces crystallization. 

Namely, Q min is given by 

_ p d c P (a + h) 3 (T c -T ) 
^ min " A;[(a + /i)3-a 3 ] " l } 

We get Q min = 0.8 x 10 27 Kcm- 3 for r coo i = 10~ 6 s and T c = 2000 K, which is the 

graphitization temperature estimated from Eq. (TT5|) . because graphitization of the carbon 

mantle leads to silicate crystallization in Exp. 2. Note that this Qmm-vahie is in good 

agreement with the numerical calculations given in section 13.11 



4. Conditions of Low Temperature Crystallization in Astrophysical 

Environments 

The previous section focused on the analysis of Exp. 2, in which partial crystallization 
of a silicate core was observed when the particle was exposed to the ambient gas of its 
pressure of latm (i.e., n ~ 10 19 cm -3 ). For low gas densities as those in astrophysical 
environments, however, crystallization becomes easier than in the conditions in Exp. 2 
because of inefficient collisional cooling by the ambient gas. Radiative cooling dominates 
the collisional cooling for very low gas densities. The gas density for which radiative and 
collisional coolings become comparable is estimated from A ra( j = A co u to be 

1/2 



n 



kv V0.01/ \2.2/ V T o J ViOOOKy v ; 



for T 3> Tp, where ji = 2.2 is the mean molecular weight of a gas of the solar composition 



( lAsplund et al. 



20061 ). As will be shown later (Figs. [7]and[H]), complete crystallization of 
the silicate core is realized for the gas densities lower than that estimated above, unless 
the stored energy density Q is so small or the reaction timescale r is so long to realize the 
crystallization. 
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We consider two possible processes of crystallization induced (A) by graphitization 
followed by chemical reactions in the mantle, and (B) reactions alone. The case (B) is 
considered to clarify the effect of the reactions without pre-heating played by the amorphous 
carbon mantle in Exp. 2. Actual situations may be situated between these two possibilities 
depending on the amorphousness and the composition of the mantle. 

First we consider the case (A). We carry out numerical calculations with varying the 
ambient gas density, the stored energy density Q, and the reaction timescale r and explore 
the crystallization conditions. Figure [7J summarizes the result as a function of the number 
density of ambient gas molecules, showing the ranges of Q (the upper panel) and r (the 
lower panel), in which amorphous silicate particles crystallize. It should be noted that, as 
the gas density decreases, the regions of Q and r inducing the crystallization widely extend 
to small Q and long r. This extension is caused by the slow cooling due to inefficient 
collision frequency of the ambient gas. Complete crystallization of the silicate core is 
possible for the gas density n < 10 18 cm -3 (depending on the stored energy density Q and 
the timescale of the reactions r); the gas density in Exp. 2 is found to be marginal and 
results in partial crystallization of the silicate core. 

Next we explore the case (B) of crystallization induced by chemical reactions alone. 
Figure M shows the result of the numerical calculations. The overall tendency is the 
same as in the case (A). Namely, the ranges of Q and r yielding the crystallization 
extend to the regions of small Q and long r with decreasing the gas density. Compared 
to the case (A) shown in Fig. [7J quantitatively, the minimum value of Q required for 
the crystallization increases about twice of that in the case (A). In addition, partial 
crystallization occurs even for low gas densities. For n < 10 13 cm -3 , partial crystallization 
occurs in the range of 0.8 < Q < 1.1 x 10 27 Kcm -3 , while complete crystallization occurs 
for Q > 1.1 x 10 27 Kcm" 3 . 
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We formulate the crystallization conditions that reproduce the numerical results shown 
in Figs. [7] and El The condition for the crystallization to occur is given by 



Q > Qmm and r < r cool , 



(25) 



where r coo i is the cooling timescale of the grain and Q m i n is the minimum stored energy 
density to induce crystallization. The minimum stored energy density Q m i n is given by 
Eq. (|23|) . in which T c depends on the case (A) or (B). In the case (A), crystallization of the 
amorphous silicate core is induced by graphitization of the amorphous carbon mantle, and 
the temperature T c is estimated from the graphitization condition (fl~5]) expressed as 

-l 



rp _ E c 
J- r 



In 



2fic /3 z/ c p d c p 



? + l 



(T c » T ). 



3(knv + 2a B e c T^) \h 

The solid line in the upper panel of Fig. [7] shows Q mm with T c determined from Eq. ( J26l) 
In the case (B), T c is estimated from the relation that 

2f\ 1 i /3 z/ si p d c p 



(26) 



k 



In 



I - 
3(knv + 2a-QE c T^) \ a 



(27) 



derived from the silicate crystallization condition that 



7~cool • 



(28) 



\da B i/dt\ 

The solid line in the upper panel of Fig. [S] shows Q m i n using Eq. (j27|) . One should note 
that both Q m in thus estimated well reproduces the results of the numerical calculations. It 
should be pointed out that Q m i Q does not depend on the absolute size of a grain but on the 
ratio between the mantle thickness h and the core radius a. The value of Q m ; n decreases 
with increasing the ratio h/a. The cosmic abundance o f elements puts a constraint on the 



ratio h/a to range from 0.3 to 0.5 (JKimura et al 



20031 ). This implies that Q m \ n does not 



depend much on h/a in astrophysical situations. We used h/a = 0.4 as a standard value in 
this study. 
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A maximum value of reaction timescale r allowing crystallization is given by r coo i (see 
Eq. (I25p). The dot- dashed lines in the lower panels of Figs. [7] and [S] show r coo i given by 
Eq. ( 1T2|) as a function of the gas density n. The cooling time becomes constant for n < 10 13 
cm -3 , for which density the cooling is almost due to thermal radiation from the particles. 
The difference in r coo i is small in cases (A) and (B). One should note that r coo i reproduces 
the boundary of the crystallization region in the n-r plane, indicating that the condition 
Eq. (J25J works. 

5. Discussion 

We have constructed a theoretical model for low-temperature crystallization induced 
by exothermic chemical reactions in submicrometer-sized grains having a silicate core and a 
carbonaceous mantle. The validity of our model is proved by its application to the results 
of laboratory experiments on crystallization of amorphous silicate grains by Kaito and his 
colleagues. This enables us to formulate crystallization conditions of amorphous silicate in 
low-temperature environments using the stored-energy density Q and the reaction timescale 
r, irrespective of details of chemical reactions (see (1251) ). It should be pointed out that the 
conditions given in Eq. (12"5|) are valid regardless of the grain structure, in spite of the fact 
that the equations are derived from the assumption of a core-mantle structure. We conclude 
that the low-temperature crystallization mechanism discussed in this paper should work at 
any dust grains unless the density Q of energy stored in the grains lies below the critical 
density Q min . 

Our results suggest that crystallization of amorphous silicate takes place in much 
lower temperatures than hitherto considered. As stated in section 1, crystallization of 
amorphous silicate needs high temperatures ( t 1000 K), while the present mechanism 
works by moderate heating, say, of a few hundred kelvin, depending the gas density of 
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interest. In a protoplanetary disk, the number density of gas and i ts temperatur e at 1 



AU is estimated to be n ~ 10 14 cm 3 and T = 300 K, respectively (JHayashi 



19811). One 



could derive the conditions for inducing silicate crystallization from Figs. [7] and [H] to be 
Qmin = (4 — 8) x 10 26 Kcm~ 3 and r coo i = 10~ 2 s. This Qmin-value corresponds to the 
concentration of reactive molecules in the mantle being (5-10) % for the heat of chemical 
reaction E T ~ 10 5 K. This concentration of reactive molecules is in harmony wi th (1-10) % 



suggested in a carbonaceous mantle of interstellar dust grains ( IGreenberg 



19761 ). Therefore, 



the nonthermal crystallization of silicate grains might come into effect in protoplanetary 
disks if heat of chemical reactions is deposited on the grains within a short timescale. 

The minimum stored energy Q mm measures hardness of crystallization of silicate of 
various composition in the nonthermal crystallization; the larger Q mm , the harder the 
crystallization is. When the crystallization is induced by reactions (case (B) in section 
4), Qmin is determined mainly by the activation energy of crystallization of the silicate 
(see Eq. (f23|) and Eq. fTSTl) ). For example, the activation energies of crystallization of the 
amorphous silicate are E^ = 39000 K for 



enstatite composition (IFabian et al. 



brsterite composition and E s i = 42000 K for that of 
200oJ). The Q min values are Q min = 8 x 10 26 Kcm- 3 for 
forsterite composition and Q m [ n = 9 x 10 26 Kcm -3 for enstatite composition. This indicates 
that the crystallization of amorphous silicate of enstatite composition is harder than that 
of forsterite composition. On the other hand, in the case that silicate crystallization is 
induced by crystallization of the amorphous mantle (graphitization for an amorphous 
carbon mantle) due to deposition of heat of reactions in the mantle (case (A)), Q m { n is 
determined by the activation energy of crystallization of the mantle material. In this case, 
the hardness of silicate crystallization is independent of its composition. 

As long as there is a mechanism to trigger chemical reactions in a carbonaceous 
layer of dust grains, amorphous silicate in the grains would suffer from low-temperature 
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crystallization. One of the plausible mechanisms is a flash heating of dust by shock waves, 
which has been proposed as a formation mechanism for chondrules found in meteorites. 
Once the dust is heated moderately in a shock region, diffusion of reactive molecules 
becomes active and will lead to chain chemical reactions. Note that the heating required 
for inducing crystallization of amorphous silicate is much lower than the one for chondrule 
formation. Collisions among dust particles, which raise the temperature in a local region of 
their surfaces (~ tens of degree), may be another mechanism of flash heating. Since the 
above-mentioned mechanisms are most likely common to happen in protoplanetary disks, 
such a flash heating would trigger silicate crystallization in a protoplanetary disk. 

Even if shock waves and mutual collisions are not effective enough to trigger chemical 
reactions, the low-temperature crystallization of amorphous silicate could be induced by 
spontaneous heating k nown as the Wigner storage mechanism in nuclear reactor engineering 



(jShabalin et al. 



20031 ). It is well known for graphite irradiated by neutrons that energy 



stored in the form of lattice defects will be released spontaneously when the number density 
of the defects reaches a certain value. A similar phenomenon of spontan eous stored-ener gy 



release was observed i n experimenta 



Shabalin et al. 



2003 



Kulagin et al. 



stud ie s of neutron-irra diated ices (ICarpenter 



1987 



20041 ) . ICarpenterl ( 119871 ) found that the temperature 



of a cold solid methane moderator rises rapidly and uncontrollably when the heat loss 
becomes less effective after irradiation of fast neutrons into the moderator. The rapid 
temperature rise is due to the phenomenon that the energy stored at low temperatures 
in the form of reactive molecules was abruptly released from the irradiated methane by 
thermally activated diffusion of the molecules and their subsequent reactions. 

Reactive molecules in the ice mantle of grains would form by irradiation of UV or high 
energy particles such as cosmic ray in the outer cold region of a protoplanetary disk or in 
a molecular cloud. When reactive molecules in such a grain are accumulated to a critical 
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number density, the heat released by the activation of reactions would induce crystallization 
of amorphous silicate in the grain. However, there is no observational evidence for the 
presence of crystalline silicate in the outer cold region of a protoplanetary disk nor in a 
molecular cloud. This implies either that there is an observational selection effect to obscure 
the presence of crystalline silicate in those environments or that the temperatures are too 
low to activate chain chemical reactions in the outer cold region of a protoplanetary disk 
and in a molecular cloud. If the latter is the case, observations with high spatial resolutions 
would constrain the conditions to trigger silicate crystallization and the accumulation rate 
of reactive molecules in low-temperature astrophysical environments. 

We are grateful to C. Kaito and Y. Kimura for helpful comments and discussion and an 
anonymous reviewer for useful comments to improve the paper. T. Y. acknowledges support 
by the Grant-in- Aid for Scientific Research (21244011) from JSPS. H. K. is supported by 
the grants from CPS, JSPS, and MEXT Japan. 

A. Determination of the Activation Energy of Crystallization 

Crystallization of an amorphous solid occurs through nucleation and growth of a 
crystalline nucleus. Empirically, a characteristic time scale r of crystallization of amorphous 
silicate placed at temperature T is defined by 

r = tf«p(;§), (Al) 



where E is the effective activation energy (IHallenbeck fc Nuth 



1998 



Fabian et al. 



20001 ). 



There are two kinds of nucleation: one is homogeneous nucleation and the other is 
heterogeneous nucleation in which impurities or contact interfaces act as seed nuclei for 
crystal growth. Homogeneous nucleation occurs in the absence of the impurities or the 
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interfaces. The timescale of crystallization due to homogeneous nucleation is controlled by 
the ra te of nucleation of c rystalline nuclei, J, and the growth rate w of the crystal, and is 



given (JKouchi et al. 



19941 ) by 



Tho 



1/4 



irJw 3 
where J is a function of the interfacial energy a between amorp hous solid and crystals 



(A2) 



and the actiyation energy of mo 



Kouchi et al. 



1994 : 



ecular diffusion E si is expressed (ISeki fc Hasegaawa 



Tanaka et al. 



1981 



2008|) as 



J 



47TZA 



(3a 



3fi s ; V vr/cT 



E s j 
exp < —777; 



A{l3afT 



(A3) 



kT 27ql si k(AT si ) 2 

Here, AT si = T c si — T is the supercooling, (3 is the geometrical parameter depending on the 
shape of crystalline nuclei (4.8 for spherical nuclei). The growth rate of the crystal w is 
given by Eq. ([7j), i.e., w = —da s \/dt. 

For heterogeneous nucleation, the timescale of crystallization depends on the number 
of the nuclei in the amorphous solid. A picture of amorphous solid is that it consists of 
an ensemble of microcrystals, which have sizes of the lattice scale and are distributed in 
random orientations. According this view, there are already enough number of crystalline 
nuclei, and the timescale of (macroscopic) crystallization is given by 

ft 1 / 3 



7~hctcro 



W 



(A4) 



In this case, E s \ corresponds to the effective activation energy E in Eq. ( 1A1I) if one ignores 
the melting ( e ~ qi > siATsi/kT2 < 1). 

Figure |9] shows the crystallization timescales as a function of temperature. From 
the figure, we put constraints on the values of E si and a . Namely, we find that 
E si /k ~ 38000 — 39000 K and a < 400 erg cm -2 satisfy the results of the experiments for 
homogeneous nucleation, while E si /k ~ 38500 K satisfies the results of the experiments for 



-25- 



heterogeneous nucleation. Note that the close value of the activation energy is yielded if we 
assume homogeneous nucleation or heterogeneous nucleation. This is because the timescale 
of the crystallization is mainly controlled by crystalline growth rather than nucleation in 
the relevant temperatures. The present value of E S i is i n good agreement w ith the effective 



activation energy (E s i/k = 39100 ± 400 K) estimated by 



Fabian et al 



(boom . 
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Fig. 1. — Schematic illustration of the crystallization experiments of particles having amor- 
phous silicate core and amorphous carbon mantle. The particles used in Exp. 1 are those 
consisting of an amorphous silicate of forsterite composition and an amorphous carbon layer. 
Particles used in Exp. 2 are t he same as those i n Exp. 1 but contain reactive molecules 



(methane in the experiment by 



Kaito et al. 



(120071 1) in the mantle. 
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Fig. 2. — Comparison of sublimation timescale r su b (dotted-dashed line) and cooling 
timescale r coo i (solid line) of a grain placed at the ambient gas pressure of 1 atm and in 



vacuum. 
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Fig. 3. — Time variations of temperature (upper panel) and thicknesses of silicate crystal Z s i 
and graphite l c (lower panel) for the setup of Exp. 1. The ambient temperature is T = 870 K 
and the activation energy of graphitization is taken to be E c /k = 23000 K. The small panel 
(b) is an enlargement of the part of the time interval of 0.076 < t < 0.08 s. 
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Fig. 4. — Time variations of the temperature and the thicknesses of silicate crystal and 
graphite for the setup of Exp. 1 conditions with varying values of the activation energy of 
graphitization, E c , and the ambient temperature T . 
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Fig. 5. — Time variations of the temperature and the thicknesses of silicate crystal / s ; and 
graphite l c for the setup of Exp. 2. In this figure we set r = 5 x 10~ 8 s and varied the 
Q-values: (a) Q = 1.1 x 10 27 Kcm- 3 , (b) 1.0 x 10 27 Kcm" 3 , and (c) 0.9 x 10 27 Kcm- 3 . The 
bars with arrows indicate the range of the thickness of the crystalline silicate layer observed 
in Exp. 2 (Kaito et al. 2007). 
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Fig. 6. — The same as in Fig. 5 but for Q = 1.0 x 10 27 Kcm 3 and (a) r = 5 x 10 9 s, 
(b) 5 x 10~ 8 s, and (c) 2 x 10~ 7 s. 
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Fig. 7. — Results of calculations for nonthermal crystallization. Filled circles show complete 
crystallizations, triangles show partial crystallizations, and marks x show no crystallization 
of amorphous silicate. 



-36- 



E 
o 

o 



O 



1.5 






• 






• 


AA 7 - 






• • • 

crystallization 
• • • 




• 


'x - 








• 






•/ 


x - 








• 






/x 


x - 


1 






▲ 






X 


x - 






partial crystal 




-"x 


X 


x - 








X 
X 


X 
X 
X 


X X 

X X 
no crystallization 

X X 


XXX 

1 1 1 




- X 

- X 


X 

X 


0.5 


- X 


X 

■ ■i 


X 


X 


X 


X 


x - 



10 



14 



10 



16 



10' 



10 



20 



10' 



.CO. 



10" 



number density [cm ] 



! X 


X X 


X 


X 


X 


x 1 


: X 


X X 


X 


X 


X 


x 1 


r 


:^x 


no 
X 


crystallization 
X X 


x 1 


r 


s 

• • 

• • 

• • 


X 


X 
s 

X 

s 

• 


X 
X 

N 


XXX 


r 


crystallization 
• • 




• 


• 




r 


• • 

• • 




• • A ] 
partial crystal 

• • X aa; 



10 



14 



10 



16 



10 



18 



10 



20 



number density [cm ] 



Fig. 8.— The same as Fig. [7] but for no graphitization of amorphous carbon. 



-37- 



1(T 



_cd 

O 
CO 

CD 

|10 ( 



10 



-5 

500 




Fabian et al. (2000) 
- Kamitsuji et al. (2005) 



E si /k=38500K 
Forsterite 



1000 



1500 2000 



T[K] 



t 
Tmelt 



Fig. 9. — Crystallizati on time as a f unctio n of 



of the experiments by 



Fabian et al. 



temperature T. Also p lotted are the results 



( J2000h and 



Kamitsuji et al. 



teoofih 



